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ON NON-EXISTENCE OF A ONE FACTOR INTEREST RATE 
MODEL FOR VOLATILITY AVERAGED GENERALIZED 
FONG VASICEK TERM STRUCTURES 

BEATA STEHLIKOVA and DANIEL SEVCOVICi 

Abstract. We study the generalized Fong-Vasicek two-factor interest rate model with stochas- 
tic volatility. In this model dispersion is assumed to follow a non-negative process with volatility 
proportional to the square root of dispersion, while the drift is assumed to be a general function. We 
consider averaged bond prices with respect to the limiting distribution of stochastic dispersion. The 
averaged bond prices depend on time and current level of the short rate like it is the case in many 
popular one-factor interest rate model including, in particular, the Vasicek and Cox-IngersoU-Ross 
model. However, as a main result of this paper we show that there is no such one-factor model 
yielding the same bond prices as the averaged values described above. 
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1. Overview of term structure modeling. The purpose of this paper is to 
study the generahzed Fong-Vasicek two-factor interest rate model with stochastic 
volatility. In this model the dispersion of the stochastic short rate (square of volatility) 
is assumed to be stochastic as well and it follows a non-negative process with volatility 
proportional to the square root of dispersion. The drift of the stochastic process for 
the dispersion is assumed to be in a rather general form including, in particular, 
linear function having one root (yielding the original Fong-Vasicek model, cf. [6 ) or 
a cubic like function having three roots (yielding a generalized Fong-Vasicek model 
for description of the volatility clustering, see e.g. [H]). We consider averaged bond 
prices with respect to the limiting distribution of stochastic dispersion. The averaged 
bond prices depend on time and current level of the short rate like it is the case 
in many popular one-factor interest rate model including in particular the Vasicek 
and Cox-IngersoU-Ross model. However, as a main result of this paper we show that 
there is no such one- factor model yielding the same bond prices as the averaged values 
described above. 

The term structure of bond prices (or yields) is a function of time to maturity, 
state variables like e.g. instantaneous interest rate as well as several model parameters. 
It describes a functional dependence between the time to maturity of a discount bond 
and its present price. The yield of bonds, as a function of maturity, forms a term 
structure of interest rates. Figure 11.11 shows the different shapes of term structures 
observed on the market based on data by Bloomberg. We can observe various type 
of functional dependences including non-monotone term structures (Brazil) or term 
structures having two local maxima (UK). Interest rate models are often formulated 
in terms of stochastic differential equations (SDK) for the instantaneous interest rate 
(or short rate) as well as equations for other relevant quantities like e.g. volatility of 
the short rate process. In one-factor models there is a single stochastic differential 
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Fig. 1.1. Examples of yield curves of governmental bonds: Australia, Brazil, Japan, United 
Kingdom (27th May 2008). Source: \http://www. bloomberg. com\ 



equation for the short rate. The volatihty of the short rate process is given in a 
deterministic way. They can be written in the form 

(1.1) dr = ^{t,r)dt + a{t,r)dw, 

where w is a Wiener process. We recall that a stochastic process {w{t),t > 0} 
is called a Wiener process if u;(0) = 0, every increment w{t + At) — w{t) has the 
normal distribution N{0,At), the increments w{tn) — w{tn-i), — w{t„~2), 

. . ., w(t2) — w(ti) for < ti < . . . < t„ are independent and paths of the process are 
continuous (see e.g. p]). The function /i in (|l.ip determines the trend in evolution 
of the short rate, function a the nature of stochastic fluctuations. The price of a 
discount bond P{t, r) at time t when the value of short rate is r, is known to be a 
solution of the partial differential equation 

(1.2) — + (M(i, r) Xait, r))— + - rP = 0, 

with the terminal condition P{T,r) — 1 for any r > where T > is a maturity 
of the bond. Here A stands for the so-called market price of risk. The above linear 
parabolic PDE is derived by constructing a riskless portfolio of bonds and using Ito's 
lemma for evaluating differentials of a stochastic portfolio that is balanced by buying 
or selling bonds with different maturities. We refer the reader for a comprehensive 
overview of term structure modeling to the book by Kwok 9 , Chapter 7] for details. 
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We also remind ourselves that bond prices determine interest rates r) by the 
formula P — e~^(-^~*\ i.e. 

^(^,^) = - prolog 

One of the first models of the class (jl.ip has been proposed by Oldfich Vasicek in 
[l4] . In this model, the short rate process is assumed to follow a stochastic differential 
equation: 

(1.3) dr = k{0 -r)dt + adw, 

where K,0,a > are positive constants. Here a > stands for volatility of random 
fluctuations of the short rate process. Deterministic part of the process k{9 — r) 
represents a mean reversion process with a limit 9, referred to as long term interest 
rate. The speed of reversion is given by the parameter k > 0. In this model, for a 
constant market price of risk A, the corresponding PDE for bond prices 

dP - dP fj^ d'^P 

(1-4) ar + w''-)-^)*+Tis^-^'' = » 

has an explicit solution P{t,r) satisfying the terminal condition P{T,r) = 1 for any 
r > 0. It has the form 

(1.5) P(^,r)=^(^)e-^(*)^ 

where the functions A and B can be expressed in a closed form (see, e.g. [13] or [S]): 

B{t) = ^-^-^^^''\ i^^(^) = _ (J, _ _ ^ _ Zl) _ ^B{tf. 

There is a rich variety of several other models in which the SDE for the short rate 
is given by a general process of the form: 

(1.6) dr = (a + hr)dt + ar^dw. 

This class of short rate models includes the well-known Cox-IngersoU-Ross model [5] 
with 7 = 1/2. A thorough comparison of these models is a topic of the paper by 
Chan, Karolyi, Longstaff and Sanders [3]. Using generalized method of moments they 
estimated the model (|1.6|) and they studied restrictions on parameters imposed in 
this models. Their result that the optimal value of the parameter 7 is approximately 
3/2 (which is more than previous models assumed), started a broad discussion on 
the correct form of volatility. Let us note that their result is not universal, e.g. in 
[1], using the same estimation methodology but for LIBOR rates, 7 was estimated 
to be less than unity (which means that volatility is less than proportional to short 
rate, unlike in the result due to Chan, Karolyi, Longstaff and Sanders). Approximate 
formulae for bond prices when the short rate follows (|1.6p has been developed recently 
by [4] and [10]. 

In one-factor models, term structure of interest rates is a function of the short 
rate and model parameters. However, it means that as soon as the parameters of the 
model are chosen, the term structure corresponding to a given short rate is uniquely 
determined. This is a simplification from the reality, as it can be seen in Fig. 11.21 
showing the examples from EURIBOR data. To capture this feature, two-factor 
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models are introduced. In the two-factor models there are two sources of uncertainty 
yielding different term structures for the same short rate. They may depend on the 
value of the other factor. Moreover, two-factor models have more variety of possible 
shapes of term structures. 

A general two-factor model with the factors x, y is given by the system of SDEs: 

dx — fixdt + Uxdwi, 
dy = iJLydt + aydw2, 

where correlation between dwi and d'W2 is a assumed to be constant p G [—1, 1], i.e. 
E{dwidw2) = pdt. The short rate is a function of these two factors, i.e. r — r{x, y). 

Let us denote by P{t,x,y) the price of a zero coupon bond with maturity T, at 
the time t when the values of the factors are x and y. the PDE satisfied by the bond 
price, which reads as (cf. [H Chapter 7]) 

dP dP dP 

— + - K(Jx)^ + (My - ^y^y)^ 

, , ald'^P ald^P d^P , , 

The parameters A^, , \y stand for market prices of risk corresponding to factors x and 
y, resp. The equation for the bond price is equipped by the terminal condition at 
t = T, P{T, x,y) = 1 for any x, y. 

There are several ways of incorporating the second stochastic factor. Based on 
experimental data from real financial market it is reasonable to make an assumption 
that the market changes the volatility of the underlying process for the short rate. 
More precisely, the volatility of the stochastic process for the short rate is stochastic 
as well. An empirical confirmation of such an assumption can be found e.g. in the 
recent paper by the authors |13j. In the so-called two- factor models with a stochas- 
tic volatility we allow the volatility to have a stochastic behavior driven by another 
stochastic differential equation. As an example of such a two-factors model one can 
consider the Fong-Vasicek model (cf. [B]) in which the stochastic volatility follows 
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a mean reverting Bessel-square root process. Another possibility is to consider the 
generalized Fong-Vasicek model in which the drift function is no longer linear but it 
may have cubic like behavior having three distinct roots representing possible steady 
states for dispersion. By this we can model so-called volatility clustering phenomenon 
observed in real markets (see [13] for details). Now, as a consequence of the multi- 
dimensional Ito's lemma the corresponding equation for the bond price is a linear 
parabolic equation in two space dimensions. These spatial dimensions correspond to 
the short rate and volatility. 

Let us consider the process a Bessel square root process with a general drift 
function a, 

(1.8) dy — a{y)dt + uj^/ydw. 

It is well known that the density distribution of a stochastic process is a solution to 
the Focker-Planck partial differential equation (see [S]). Recall that the cumulative 
distribution function F = F{y,t) = Prob{y{t) < y\y{0) = yo) of the process y = y{t) 
satisfying (|1.8p and starting almost surely from the initial datum yo can be obtained 
from a solution F — OF /dy to the so-called Focker-Planck equation for the density 
function: 

(1-9) §? = T^(^^~^"^("(^)^^)' /(2^'0) = %-yo). 

Here 5{y — yo) denotes the Dirac delta function located at yo. The limiting density 
f{y) = hmt_^oo fiUit) of the process is therefore a stationary solution to the Focker- 
Planck equation (II. 9|) . A stationary solution satisfying f{y) = for y < is therefore 
a solution to the differential equation (see [T^) 

(1-10) Y^(j//)-a(y)/ = 0. 

Concerning structural assumption made on the drift function a : i? — > i? we shall 
henceforth assume the following hypothesis: 

(A) a is a function on [0, oo), ^^^^ > 1, limsup^^^ < 0. 

V y-*co y 

Now it follows from the assumption (A) made on the drift function a and [13l 
Lemma 2] that the stationary Focker-Planck equation (jl.lOp has a solution / that 
can be explicitly expressed as: 

m ^ Cy~^ exp (-^ --fd^) = Cy'-^-^ exp (-^ f ai^) 

for y > and f{y) ~ for y < 0. Here d(y) = {a{y) — a{Oj)/y and C > is a 
normalization constant such that f(y)dy ~ 1. For example, if we consider a mean 
reverting Bessel square root process for the stochastic dispersion y, i.e. the drift 
function is linear a(y) = Ky(0y — y) then the limiting density distribution function 
/ is the Gamma distribution function with shape parameters ^KydyjuP' and 2Ky/uj^ 
(see Kwok 0). 
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2. Averaging with respect to stochastic volatiHty and relation to one- 
factor models. Knowing the density distribution / of the stochastic volatility we 
are able to perform averaging of the bond price and the term structure with respect to 
volatility. Unlike the short rate which is known from the market data on daily basis, 
the volatility of the short rate process is unknown. The exact value of the stochastic 
volatility is not observable on the market, we can just observe its statistical properties. 
Therefore such a volatility averaging is of special importance for practitioners. 

We shall consider the following model with stochastic volatility: 

(2.1) dr ^ K{e - r)dt + y/ydwr, 

(2.2) dy — a{y)dt + ujy^dwy, 

with uncorrelated increments dwr and dwy of a Wiener process. The market prices 
of risk are assumed to have a form \r{r,y) ~ and \y = ^^/y- Then the bond 
price Tr{t,r,y) satisfies the following PDE: 

(2.3) - + - r) - A,)- + [a[y) -Xy)- + -y— + - r. ^ 

with the terminal condition 7r(T, r,y) = 1 for any r, y > 0. The explicit solution can 
be written as 

(2.4) ^(i,r,y)=A(t,y)e-^M'- 

with the terminal conditions A{y,T) — 1 for any y > and B{T) — 0. The solution 
can be obtained by solving the following differential equations for the functions A and 
B: 

(2.5) -B' + kB -1 = 0, 

, , dA , y ^ ^ , , ^ ~dA d'^A 

(2.6) — - B{nd - Ay - \b)A + (a(y) - Ay)^ + -y^ = 0. 

In what follows, we shall denote by (■0) the averaged value of the function t/j : 
[0, cxd) — + R with respect to the limiting density /, i.e. (?/;) = 4'{y)f{y) dy, where 
/(y) satisfies the stationary Focker-Planck equation (|1.10p The averaged bond price 
with respect to the limiting distribution of the volatility is given by 

(2.7) P{t,r) = {^{t,r,.))=a{t)e-B^'>, 

where a{t) = A(t,y)f{y)dy. 

The function P = P{t, r) is a function of time t and short rate r. Notice that it is 
the same functional dependence as for the bond price in one factor models, including, 
in particular, a solution to (|1.3p given by (jl.Sp . However, we show that there is no 
such one-factor model yielding the same bond prices as those of averaged bond prices 
P{t, r) given by the averaging of the two-factor model. 

Now we are in a position to state our main result of the paper. We are going 
to prove that there is no one-factor interest rate model for the corresponding bond 
prices that is identical with the volatility averaged bond price P{t, r) for any t £ [0, T] 
and r > 0. Suppose to the contrary that P = P{t, r) is a bond price from one-factor 
model in which the short rate is assumed to follow a general SDE 



(2.8) 



dr = fi(r)dt + Q{r)dw. 
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Then is satisfies the PDE 

BP BP ^ f)^P 

(2.9) — + (M(r) - A(r))— + ^^^r)— -rP = 0, 

where ^(r) is the drift of the short rate process, ri(r) is its volatility and A(r) is the 
product of the volatility and the corresponding market price of risk of the model. 
Substituting the form of the solution we obtain that 

(2.10) + 

We see that the left hand side is a function of t only. We denote it by (j){t), i.e. 
(I)(t) = a{t)/{a{t)B{t)). Then, 

(2.11) nr + ^i{r) - A(r) = 0(<) + ^n'^{r)B{t), 



and so the right hand side is constant with respect to t. Hence for any t it equals 
(j){T) + if22(r)B(r) = (/)(r), which is a constant denoted by K. We have 

(2.12) Kr + fi{r) ~ A(r) = ^n\r)B{t) + (j){t) = K, 

from which it follows that fl{r) = f2 is a constant and that /i(r) — A(r) ^ K — nr. 

Let us denote by and d the first two statistical moments of the random variable 
with respect to the limiting density function /, i.e. 

/•oo poo 

= {y)= / yf{y)dy, d = {y') = / y^f{y)dy. 
Jo Jo 

We know that a{T) — 1 and from the expression 

(2.13) a'it) = (k- ^B{t)] a{t)B{t) 



2 

we can recursively compute the values of the time derivatives of function a at time T: 

(2.14) a'(T) = 0, 

(2.15) a"(T) = -K, 

(2.16) a"'{T) = -Kn-n^, 

(2.17) a""{T) = 3K^ - SCl^n ~ Kn. 

Another way of computing these derivatives is using the expression 

/•oo 

a{t)^ / A{t,y)f{y)dy 



and the partial differential equation (|2.6p for the function A and the stationary Focker- 
Planck equation for the limiting density function /. Indeed, by integration by parts 
and taking into account boundary conditions yf{y) = for y = 0, +oo for the limiting 
density / we obtain 



I{y)dy 



.OA , uj'^ d'^A " 
dy'^ 

r°° ( ri I P- rP- \ 



p / , .dA ^2 ^ 
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Furthermore, by p.lOp . we have 

f°° dA 2 f°° 

/ -T-yf{y)dy = — 2 / Aoii.v)f{y)dv- 

Jo cly ^ Jo 

Therefore 

a'{t) = (^{t) (ne -Xy- - ^c^iy)^ A{t, y)f{y)dy. 

Now taking into account PDE (|2.6p for the function A{t, y) we can recurrently evaluate 

OA d'^A 
A{T,y) = l, _(T,y) = 0, —{T,y) = -Ke + Xy, 

d'^A 

-^{T, y) = -n^e - (1 - KX)y - X{a{y) - Xy), 

for any y > 0. Using the above expressions and the identities B{T) — 0, B'{T) = 
— 1,B"{T) = —K,B'"{T) = — K^, after straightforward computations we obtain 

(2.18) a'(T) = 0, 

(2.19) a"{T) = -k9 + Xa^ 

and by comparing (|2.15p and (|2.19p we obtain the expression for the constant K in 
terms of the model parameters as 

(2.20) K = K9-Xa'^. 
Computing the next derivative we end up with 

(2.21) a"'{T) = XXa^ -K^9 + KXa^ -cr^ 

and by comparing (|2.16p and (|2.2ip we can express the volatility ft as 

(2.22) ^a^{l-XX). 

Notice that the PDE for the averaged bond price now reads as follows: 

9P , ,n ^ ^ 2^ a^il-XX)d^P „ 

which is the PDE corresponding to the classical one-factor Vasicek interest rate model. 
Now we fully determined the drift function minus market price function fi{r) — A(r) 
as well as the volatility function fl{r) in (|2.9p . 

In order to achieve contradiction we finally compute the fourth derivative as 

a""{T) = 3X^d + {-GkOX + k^X - 3k + X{kX - 1 + AA))^^ 



2 r. 



(2.23) +3K^r - + ^AA / a\y)f{y)dy. 

^ Jo 

Comparing (|2.17p and (|2.23p we get the condition 

(2.24) cr2(2KAA + 1 - AA^) = — AA / a^{y)f(y)dy + 3X^(d-<T'^) 

^ Jo 



On non-existonce of a one factor interest rate model for volatility averaged term structures 9 



However, the latter equality can not be satisfied for general choice of model parame- 
ters. Indeed, setting A = 0, we obtain = = yf{y)dy which is not possible as 
f{y) > for y > 0. 

Summarizing, we have shown the following theorem: 

Theorem 2.1. Consider the generalized Fong-Vasicek two-factors model with 
stochastic volatility \2. 2\) and the averaged bond price P{t,r) with respect to the limit- 
ing distribution of the stochastic dispersion. Then there is no one-factor interest rate 
model 112. 8\) with corresponding PDE for the bond price 112. 9\) yielding the same bond 
prices as the averaged values P{t, r) from the two-factor model. 
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